1 

Error Probability Bounds for 
Balanced Binary Relay Trees 

Zhenliang Zhang, Student Member, IEEE, Ali Pezeshki, Member, IEEE, 
William Moran, Member, IEEE, Stephen D. Howard, Member, IEEE, 
and Edwin K. P. Chong, Fellow, IEEE 



Abstract 

We study the detection error probability associated with a balanced binary relay tree, where the leaves of 
the tree correspond to N identical and independent detectors. The root of the tree represents a fusion center that 
makes the overall detection decision. Each of the other nodes in the tree are relay nodes that combine two binary 
messages to form a single output binary message. In this way, the information from the detectors is aggregated 
into the fusion center via the intermediate relay nodes. In this context, we describe the evolution of Type I and 
Type II error probabilities of the binary data as it propagates from the leaves towards the root. Tight upper and 
lower bounds for the total error probability at the fusion center as functions of N are derived. These characterize 
how fast the total error probability converges to with respect to N, even if the individual sensors have error 
probabilities that converge to 1/2. 

Index Terms 

Binary relay tree, distributed detection, decentralized detection, hypothesis testing, information fusion, dynamic 
system, invariant region, error probability, decay rate, sensor network. 

I. Introduction 

Consider a hypothesis testing problem under two scenarios: Centralized and decentralized. Under the 
centralized network scenario, all sensors send their raw measurements to the fusion center which makes 
a decision based on these measurements. In the decentralized network introduced in [1J, sensors send 
summaries of their measurements and observations to the fusion center. The fusion center then makes 
a decision. In a decentralized network, information is summarized into smaller messages. Evidently, the 
decentralized network cannot perform better than the centralized network. It gains because of its limited 
use of resources and bandwidth; through transmission of summarized information it is more practical and 
efficient. 

The decentralized network in [1J involves a parallel architecture, known as a star architecture [0Q- 
lfT5ll . ll31"ll . in which all sensors directly connect to the fusion center. A typical result is that under 
the assumption of (conditionally) independence of the sensor observations, the decay rate of the error 
probability in a parallel network is exponential (61 . 

Several different sensor topologies have been studied under the assumption of conditional independence. 
The first configuration for such a fusion network considered was the tandem network [|T6Tl — [|20Tl . fl3TTl . In 
such a network, each non-leaf node combines the information from its own sensor with the message it 
has received from the node at one level down, which is then transmitted to the node at the next level up. 
The decay rate of the error probability in this case is sub-exponential 11201 . As the number of sensors N 
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goes to infinity, the error probability is Q(e~ cN ) for some positive constant c and for all d G (1/2, 1) 
|fT8ll . This sensor network represents a situation where the length of the network is the longest possible 
among all networks with N leaf nodes. 

The asymptotic performance of single-rooted tree networks with bounded height is discussed in ETTl - 
ll29llll3T1l . Even though error probabilities in the parallel configuration decrease exponentially, in a practical 
implementation, the resources consumed in having each sensor transmit directly to the fusion center might 
be regarded as excessive. Energy consumption can be reduced by setting up a directed tree, rooted at the 
fusion center. In this tree structure, measurements are summarized by leaf sensor nodes and sent to their 
parent nodes, each of which fuses all the messages it receives with its own measurement (if any) and 
then forwards the new message to its parent node at the next level. This process takes place throughout 
the tree culminating in the fusion center, where a final decision is made. For a bounded-height tree, the 
error exponent is as good as that of the parallel configuration under certain conditions. For example, for 
a bounded-height tree network with linv^-Kx, £n/tn = 1, where tn denotes the total number of nodes 
and tjsi denotes the number of leaf nodes, the optimum error exponent is the same as that of the parallel 
configuration j|22|. 

The variation of detection performance with increasing tree height is still largely unexplored. If only the 
leaf nodes have sensors making observations, and all other nodes simply fuse the messages received and 
forward the new messages to their parents, the tree network is known as a relay tree. The balanced binary 
relay tree has been addressed in QUI . Evidently, the height of this tree is log 2 N. In ll30l . it is assumed that 
the leaf nodes are independent sensors with identical Type I error probability (a ) and identical Type II 
error probability ((3 ). It is shown there that if sensor error probabilities satisfy the condition a + (3 < 1, 
then both the Type I and Type II error probabilities at the fusion center both converge to as the iV goes 
to infinity. If a + A) > 1> then both Type I and Type II error probabilities converge to 1. 

We consider the same tree configuration in this paper and describe the precise evolution of Type I and 
Type II error probabilities in this case. In addition, we provide upper and lower bounds for the total error 
probability at the fusion center as functions of N. These characterize the decay rate of the total error 
probability. We also show that the total error probability converges to under certain condition even if 
the sensors are asymptotically crummy, that is, ao + A) — > 1- 

The rest of the paper is organized as follows. In Section n, we formulate the decentralized detection 
problem in the setting of balanced binary relay trees. In Section III, we discuss the evolution of Type 
I and Type II error probabilities. In Section IV, we derive upper and lower bounds for the total error 
probability at the fusion center as functions of N. In Section V, we discuss some corollaries focusing on 
the asymptotic regime as N — > oo. Finally, Section VI contains the concluding remarks. 

II. Problem Formulation 

We consider the problem of binary hypothesis testing between H and Hi in a balanced binary relay 
tree. Leaf nodes are sensors undertaking initial and independent detections of the same event in a scene. 
These measurements are summarized into binary messages and forwarded to nodes at the next level. Each 
non-leaf node with the exception of the root, the fusion center, is a relay node, which fuses two binary 
messages into one new binary message and forwards the new binary message to its parent node. This 
process takes place at each intermediate node culminating in the fusion center, at which the final decision 
is made based on the information received. 

In this configuration, as shown in Fig. [TJ the closest sensor to the fusion center is as far as it could 
be, in terms of the number of arcs in the path to the root. In this sense, this configuration is the worst 
case among all N sensor relay trees. Moreover, in contrast to the configuration in [|22l discussed earlier, 
in our balanced binary tree we have \im TN ^ OCJ £ N /t n = 1/2 (as opposed to 1 in [22]). Hence, the number 
of times that information is aggregated here is essentially as large as the number of measurements (cf., 
(|22l. in which the number of measurements dominates the number of fusions). 

We assume that all sensors are independent given each hypothesis, and that all sensors have identical 
Type I error probability a and identical Type II error probability (3 . The likelihood-ratio test 021 is 
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Fig. 1. A balanced binary relay tree with height fc. Circles represent sensors making measurements. Diamonds represent relay nodes which 
fuse binary messages. The rectangle at the root represents the fusion center making an overall decision. 



applied as the fusion rule at the intermediate relay nodes and at the fusion center. We are interested in 
following questions: 

• What are these Type I and Type II error probabilities as functions of iV? 

• Will they converge to at the fusion center? 

• If yes, how fast will they converge with respect to iV? 

Fusion at a single node receiving information from the two immediate child nodes where these have 
identical Type I error probabilities a and identical Type II error probabilities p. provides a detection with 
Type I and Type II error probabilities denoted by (a', /?'), and given by ll30l : 

f (l-(l-«) 2 ,/3 2 ), a</3, 
(a', 13') = f(a,/3):= I (1) 

{ (a 2 ,l - (1 -/3) 2 ), a>p. 

Evidently, as all sensors have the same error probability pair (ao,Po), all relay nodes at level 1 will 
have the same error probability pair (a x , pi) = f(a , /3 ), and by recursion, 

(a k+1 ,(3 k+1 ) = f(a k ,/3 k ), k = 0, 1, . . . , (2) 

where (a k , (3 k ) is the error probability pair of nodes at the A;th level of the tree. 

The recursive relation ([2]) allows us to consider the pair of Type I and II error probabilities as a 
discrete dynamic system. In [|30ll . which focuses on the convergence issues for the total error probability, 
convergence was proved using Lyapunov methods. The analysis of the precise evolution of the sequence 
{(a k , (3 k )} and the total error probability decay rate remain open. In this paper, we will establish upper and 
lower bounds for the total error probability and deduce the precise decay rate of the total error probability. 

To illustrate the ideas, consider first a single trajectory for the dynamic system given by equation ([T|), 
and starting at the initial state (a ,p ). This trajectory is shown in Fig. |2j It exhibits different behaviors 
depending on its distance from the (3 = a line. The trajectory approaches (3 = a very fast initially, but 
when (ak, Pk) approaches within a certain neighborhood of the line p = a, the next pair (ak+i, Pk+i) 
will appear on the other side of that line. In the next section, we will establish theorems that characterize 
the precise step-by- step behavior of the dynamic system ([2]). 
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Fig. 2. A trajectory of the sequence {(ak,Pk)} in the (a, f3) plane. 

III. Evolution of error probabilities 

The relation ([I]) is symmetric about both of the lines a + (3 = 1 and (3 = a. Thus, it suffices to 
study the evolution of the dynamic system only in the region bounded by a + (3 < 1 and (3 > a. Let 
hi := {(a, (3) > 0\a + (3 < 1 and (3 > a} be this triangular region. Similarly, define the complementary 
triangular region £ := {(a, (3) > 0\a + f3 < 1 and (3 < a}. 

LetB 1 := {(a,/3) G U\(l -a) 2 + /3 2 < 1}. If (a k ,/3 k ) G B u then the next pair (a k+1 , fi k+1 ) = f(a k ,/3 k ) 
crosses the line (3 = a to the opposite side from (a k , (3 k ). More precisely, if (a k , (3 k ) G hi, then (a k , (3 k ) G 
B\ if and only if (a k+ i, (3 k+ \) = f(a k ,(3 k ) G C. In other words, B\ is the inverse image of C under / in 
hi. The set £>i is shown in Fig. (3^a). Fig. ^b) illustrates this behavior of the trajectory for the example 
in Fig. [2} For instance, as shown in Fig. ^b), if the state is at point 1 in B-y, then it jumps to the next 
state point 2, on the other side of (3 = a. 

LetB 2 := {(a,/3) G U\{l-a) 2 + (3 2 > 1 and (l-a) 4 +/3 4 < 1}. It is easy to show that if (a k , k ) G hi, 
then (a k , (3 k ) G B 2 if and only if (a k +i, (3 k +i) = f(a k , (3 k ) G B x . In other words, B 2 is the inverse image 
of B\ in hi under /. The regions and the behavior of / is illustrated in the movement from to point 1 
in Fig. [3tb). The set B 2 is identified in Fig. [3|a), lying directly above B\. 

Now for an integer m > 1, recursively define B m to be the inverse image of B TO _i under /, denoted 
by B m . It is easy to see that B m := {(a, (3) G U\{1 - a) 2 ^ + /3 2(m_1) > 1 and (1 - a) 2 " 1 +/3 2 " 1 < 1}. 
Notice that W = |J™ =1 B m . Hence, for any (a , f3 ) G W, there exists m such that (ao,/3o) G -B m . This 
gives a complete description of how the dynamics of the system behaves in the upper triangular region 
hi. For instance, if the initial pair (a , (3 ) lies in B m , then the system evolves in the order 

B m ^ B m —\ -> >B 2 ^Bl 

Therefore, the system will enter B\ after m — 1 levels of fusion, i.e., (a m _i, /3 m _±) G -Bi. 

As the next stage, we consider the behavior of the system after it enters B\. The image of i?i under 
/, denoted by R c , is (see Fig. [3] (a)) 

i? £ := {(a, /3) G £| v 7 ! 3 ^ + v 7 /^ > 1} (3) 

We can define the reflection of B m about the line (3 = a in the similar way for all m. Similarly, we 
define the reflection of R c about the line (3 = a to be R u ■— {(a, (3) G U\y/1 — (3 + yfa > 1}. We denote 
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Fig. 3. (a) Regions Bi, E>2, and Rc in the (a,/3) plane, (b) The trajectory in Fig. [5] superimposed on (a), where solid lines represent 
boundaries of B m and dashed lines represent boundaries of R. 

the region R u U Rc by R. Below R is shown to be an invariant region in the sense that once the system 
enters R, it stays there. For example, as shown in Fig. |3jb), the system after point 1 stays inside R. 

Proposition 1: If (ak , (3k ) G R for some k , then (a fc , f3 k ) e R for all k > k . 
Proof: First we show that Bi C R u C B 1 U B 2 . 

Notice that B\, Ru, and B 1 U B 2 share the same lower boundary /3 — a. It suffices to show that the 
upper boundary of R u lies between the upper boundary of B 2 and that of B 1 (see Fig. [4]). 

First, we show that the upper boundary of Ru lies above the upper boundary of B\. We have 

1 - (1 - v^) 2 > v 7 !- (l-«) 2 
<^=^> — a > \/2a — a 2 

a 2 + a - 2ai > 0, 




Fig. 4. Upper boundaries for B\ , B2 , and Ru . 
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which holds for all a in [0, 1). Thus, B\ C Ru- 

Now we prove that the upper boundary of R u lies below that of B 2 . We have 

(l_(l_ a )4)| > 1 _( 1 _ v ^)2 




1 - (1 -a) 4 > (2v^-a) 4 

-2(y/a - l) 2 a(-c^ + c^v 7 ^ - 1) + 4:y/a(y/a - 1) + a - 2) > 0, 



which holds for all a in [0, 1) as well. Hence, Ru C B 1 U B 2 . 

Without lost of generality, we assume that (a ko , f3 ko ) G Ru- That means (a ko , f3 ko ) G B x or (a ko , /3 ko ) G 
5 2 n If (a fc() ,/3fc ) G B u then the next pair (a ko+1 , f3 ko+1 ) is in R c . If (a feo ,/3 fc() ) G B 2 n ffo, then 
(ci!fc 0+ i, Pko+i) £ Bi C Ru and (a ko+2: fik +2) G -Re- By symmetry considerations, it follows that the 
system stays inside R for all k > k . 

m 

So far we have studied the precise evolution of the sequence {(a k , f3 k )} in the (a, (5) plane. In the next 
section, we will consider the step-wise reduction in the total error probability and deduce upper and lower 
bounds for it. 

IV. Error probability bounds 

The total error probability for a node with (a k , f3 k ) is (a k + /3 k )/2 because of the equal prior assumption. 
Let L k = a k +(3 k , namely, twice the total error probability. Analysis of the total error probability will result 
from consideration of the sequence {L k }. In fact, we will derive bounds on logL^ 1 , whose growth rate 
is related to the rate of converge of L k to 0. (Throughout this paper, log stands for the binary logarithm.) 
We divide our analysis into two parts: 

I We will study the shrinkage of the total error probability as the system propagates from B m to Bi, 
II We will study the shrinkage of the total error probability after the system enters B x . 

A. Case I: Error probability analysis as the system propagates from B m to B\ 

Suppose that the initial state (a , A)) lies in B m , where m is a positive integer and m/1. From the 
previous analysis, (a m -i,/3 m -i) G B x . In this section, we study the rate of reduction of the total error 
probability as the system propagates from B m to B x . 

Proposition 2: Suppose that (a k) f3 k ) G B m , where m is a positive integer and m^l. Then, 

1 < < 2. 

Proof: If (a k , (5 k ) G B m for m/1, then 

L k+ i = l-(l-a k ) 2 + P 2 k 
Ll (a k + P k f ■ 

The following calculation establishes the lower bound of the ratio L k+ i/L\\ 

L k+1 - L 2 k = 1 - (1 - a k f + Pi - (a k + f3 k f 
= ~2a\ - 2a k (3 k + 2a k 
= 2a k (l-(a k + f3 k ))>0, 

which holds in B m . 

To show the upper bound of the ratio L k +i/ 'L\, it suffices to prove that 

L k+1 -2L\ = \-{\- a k f + 01- 2{a k + p k f 
= -3a 2 k - 4a k p k + 2a k - f k < 0. 
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The partial derivative with respect to f3 k is 

d(L k+1 - 2L\) 



d/3 k 



-2fi k - Aa k < 0, 



which is non-positive, and so it suffices to consider values on the upper boundary of B\. 

L k+l -2L\ = \-{\- a k f + Pi - 2{a k + (3 k ) 2 
= 2p 2 k - 2{a k + (3 k ) 2 < 0. 

In consequence, the claimed upper bound on the ratio L k+ i/ ' L\ holds. The reader is referred to Fig. [5] for 
a plot of values of L k+ i/L 2 k . 
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Fig. 5. Ratio L^+i/ L k in Um=2 Each line depicts the ratio versus a for a fixed /3. 



Proposition 3: Suppose that (a , f3 ) E B m , where m is a positive integer and m ^ 1. Then, for 

i k 



k = 1, . . . , m — 1 



2 K logLn 



<logL fe 1 <2 fc logL - 1 



2 k 

Proof: From Proposition 2 we have, for k — 0, . . . , m — 2, 
for some 6 [1, 2]. Then for A; = 1, . . . , m — 1, 



L k = ^Yl °y 



where aj G [1,2]. Therefore, 



log L fc 1 = - \ lo S a * J - lo S L t = ~ [J2 lo S a + 2 * lo S L o 
Since logi^ 1 > and < loga^ < 1 for each i, we have 



logL^<2' c logL 



Finally, 



logL- 1 > -k + 2 fe logLo 1 = 2 k (logLo 1 - ^ . 
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Suppose that the balanced binary relay tree has N leaf nodes. Then the height of the fusion center 
is logiV. For convenience, let = L\ ogN be (twice) the total error probability at the fusion center. 
Substituting k = log N into Proposition 3, we get the following result. 

Corollary 1: Suppose that (a , A>) E B m , where mis a positive integer and m ^ 1. If logiV < m, then 

N (logLo 1 - < logP/ < iVlogLo 1 . 

Notice that log N/N tends to zero as N tends to infinity. Hence, the decay rate of P/y with respect to 
N is exponential, with error exponent logZ-Q 1 . In other words, as (a k , f3 k ) marches from B m towards B%, 
the total error probability at each level reduces exponentially quickly. 



B. Case II: Error probability analysis when the system stays inside R 

We have derived error probability bounds up until the point where the trajectory of the system enters 
B x . In this section, we consider the total error probability redution from that point on. First we will 
establish error probability bounds for even-height trees. Then we will deduce error probability bounds for 
odd-height trees. 

1 ) Error probability bounds for even-height trees: If (ato, (3o) E B m for some m/1, then (a m _i, m -i) E 
Bi. The system afterward stays inside the invariant region R (but not necessarily inside Bi). Hence, the 
decay rate of the total error probability in the invariant region R determines the asymptotic decay rate. 
Without lost of generality, we assume that (ckq, /3q) lies in the invariant region R. In contrast to Proposition 
2, which bounds the ratio L k+ i/ L\, we will bound the ratio Lk+2/Ll associated with taking two steps. 

Proposition 4: Suppose that (a k , (5 k ) E R. Then, 

1 < < 2. 

Proof: Because of symmetry, we only have to prove the case where (a k , (3k) lies in Ru- We consider 
two cases: (a k , (3 k ) E B x and (a k , (3 k ) E B 2 H R u . 
In the first case, 



L k+2 (1 - (1 - ajk ) a ) 2 + 1 - (1 - ^ 2 



L 2 k (a k + (3 k y 

To prove the lower bound of the ratio, it suffices to show that 

L k+2 - Ll = (1 - (1 - a k ) 2 ) 2 + 1 - (1 - dlf - (a k + (3 k f 

= (a k + h- \){{a k - (3 k ) 3 + 2a k p k (a k - f3 k ) - (a k - $ k ) 2 - 2a 2 k ) > 

which holds for all (a k , (3 k ) E B\. 

To prove the upper bound of the ratio, it suffices to show that 

L k+2 - 2L\ = a\- Aa\ + 2a\ - 4a k f3 k - (3 4 k < 0. 

The partial derivative with respect to f3 k is 

d(L k+2 - 2Ll) ^ 

^ = -Aa k - A(3i < 

which is non-positive. Therefore, it suffices to consider its values on the curve (3 k = a k , on which 
L k + 2 — 2L\ is clearly non-positive. The reader can also refer to Fig. j6^a) for a plot of values of L k+2 /L\ 
in Bl 



Now we consider the second case, namely (a*, 0k) £ B 2 D Ru, which gives 



L k+2 l-(l-a k ) A + Pt 



To prove the lower bound of the ratio, it suffices to show that 

L k+2 -L 2 k = (l-{1- a k ) A ) + f3i - (a k + k f 

= -(a k +0 k - l)(a 3 k - a\0 k - 3a 2 k + a k /3 2 k + 2a k /3 k - /3 3 k - 0\ + Aa k ) > 0. 

Therefore, it suffices to show that 

<f>(a k , 0k) =a\ - a\0 k - 3a 2 k + a k 2 k + 2a k k - 0\ - 0\ + Aa k > 0. 
The partial derivative with respect to k is 

^- = -(a k - k f - 20 2 k + 2(a k - k ) < 0. 

Opk 

Thus, it is enough to consider the values on the upper boundaries s/l — k + ^fa~k = 1 and a k + k 
If a k + 0k = 1, then the inequality is trivial, and if y/1 — k + Ja~k = 1, 



L k+2 — L k = 2(4(1 - 2^a~ k )(2a k - Q^a~ k + 5) 

and the inequality holds because ak < \ in region B 2 D Ru- 

The claimed upper bound for the ratio L k+2 / L\ can be written as 

L k+2 - 2L\ = (1 - (1 - a k f) + 0t- 2(a k + k f 

= -a\ + Aa\ - 8a 2 k + 4 «fc ~ ^k0k + 0t ~ < 0. 



' --4a k + 40 3 k -40 k <O. 



The partial derivative with respect to k is 

d(L k+2 - 2L 2 
d0 k 

Again, it is sufficient to consider values on the upper boundary of B\. Therefore, 

L k+2 - 2L\ = 20 2 k - 2(a k + k ) 2 < 0. 
The reader is referred to Fig. |6lb) for a plot of values of L k+2 /L 2 k in B 2 C\ Ru- 
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Fig. 6. (a) Ratio Lk+2/Lj, in region B\. (b) Ratio Lk+2/ Lj. in region B2 n Ru- Each line depicts the ratio versus a for a fixed j3 
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Proposition 4 gives bounds on the relationship between L k and L k+2 in the invariant region R. Hence, 
in the special case of trees with even height, that is, when log N is an even integer, it is easy to bound 
Pjv in terms of L . In fact, we will bound logP^ 1 which in turn provides bounds for P N . 

Theorem 1: If (a 0: /3 ) is in the invariant region R and logiV is even, then 

(bgLo 1 - l ^0 ) j < logP^ < v^VlogLo 1 . 

Proof: If (a , /?o) £ R, then we have (a k) /3 k ) £ R for k — 0, 1, . . . , log N — 2. From Proposition 4, 
we have 

for A; = 0, 1, . . . , log N — 2 and some a k £ [1, 2]. Therefore, for k — 2, 4, ... , log JV, we have 

A/2 \ 



where a; e [1,2]. Substituting A; = logA^ we have 

'log v^/V \ / log 



n «. *r = n * ) 

i=l / \ i=l 



Hence, 



^ log a, 1 +v / iVlogL 1 . 
Notice that logZ-Q 1 > and, for each i, < logo; < 1. Thus, 

5—1 ^ ,/^\7i„„ r — 1 



logP^ 1 ^ ViVlogLo 1 . 

Finally, 

logP/ > -logv / iV + v^logLo 1 = y/N ( log^ 1 - l °' V ' 



N 

2) Error probability bounds for odd-height trees: Next we explore the case of trees with odd height, i.e., 
log N is an odd integer. Assume that (ao, (3 ) lies in the invariant region R. First, we will establish general 
bounds for odd-height trees. Then we deduce bounds for the case where there exists (a k , /3 k ) E B 2 DRu 
for some k e {0, 1, . . . , log N — 1}. 

For odd-height trees, we need to know how much the total error probability is reduced by moving up 
one level in the tree. 

Proposition 5: If (a k , (3 k ) e U, then 

and 

Lk 

Proof: The first inequality is equivalent to 

L k+1 - L\ = 1 - (1 - a k f + f3 2 k - (a k + f3 k f 
= 2a fc (l - (a k + fa)) > 0, 
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which holds for all (a k , j5 k ) G U. 

The second inequality is equivalent to 

L k+1 - L k = 1 - (1 - a k f +Pl- (a k + (3 k ) 

= {a k - (3 k ){l-{a k + (3 k )) <0, 

which holds for all (a k , j3 k ) G W. Fig. [7] gives a plot of values of L k+ i/ L\ and L k+X /L k in W. 





Fig. 7. (a) Ratio Lt+i/L^ in region W. (b) Ratio Lk+i/Lk in region W. Each line depicts the ratio versus a for a fixed /3. 



Using Propositions 4 and 5, we are about to calculate error probability bounds for odd-height trees as 
follows. 

Theorem 2: If (a , A)) £ -R an d logiV is odd, then 



log L n 1 - 



Proof: By Proposition 5, we have 



log /7V 



< logP^ 1 < V^JVlogLo 1 . 



L\ = clLq 



for some a > 1. And, by Proposition 4, the following identity holds, 
for k = 1, 3, ... , log N — 2 and some afc G [1, 2]. Hence, we can write 

\k-l)/2 

. 2 (fc+l)/2 



L k = a \ ] [ Oj L„ 



8=1 



where 1 < «i < 2 for z = 1, 2, . . . , (k — l)/2 and a > 1. Let = log N, we have 



Pn = a 



n » 

V 1=1 / 



n » 

V 1=1 ) 



L 



/2N 
> 



and so 



log 



logPjv 1 = -logo- ^ lo§ 



i=i 



+ V^logLn 1 . 
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Notice that logL 1 > and for each i, loga; > 0. Moreover, log a > 0. Hence, 



It follows by Proposition 5 that 
for some a < 1. Thus, 



logP/^v^logLo 1 . 



L\ = clLq 



(fc-l)/2 



L k = a 



n ai i l o 



2 (fc-l)/2 



i=l 



where 1 < at < 2 for z = 1, 2 . . . , (A; — l)/2 and a < 1. Hence, 



P N = a 



n - 



V i=1 / 



L n = a 



n * 



V 4=1 / 



U 2 



and so 



/ log 



log P^ 1 = - log a 



\ 



log a* 



V 



i=l 



/ 



logLo 1 . 



Notice that logL 1 > and for each i, < logaj < 1 and log a < 0. Thus, 



log P/ > - log J j + J y log Lo 1 



2l l0gL ° 



Now we consider the special case where there exists k 6 {0, 1, . . . , log AT — 1} such that (a k , G 

P 2 n R u . 

Proposition 6: Suppose that (a k , h) £ Pi and (aifc_i, fik-x) £ P2 nP^. Then, 

I<^±1<1. 
2 " L fc " 

Proof: The right inequality is trivial. By Proposition 2, if e P 2 D Ru, then 

1 < — < 2 

-LI 



J k-1 



i.e., 



1 7~ 2 

-<%i<l, 

2 Lfc 

and in consequence of Proposition 4, if (a k -i, Pk-i) ^ B 2 f] Ru, then 

1 < ^*ti < 2 . 



Therefore, 



1 ^ Lfc+i 



2 Lfc 

In this case, Fig. M gives a plot of the ratio L k+1 /L k in the region f(B 2 U P^). 
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Fig. 8. Ratio Lk + i/Lk in the region /(£>2 U Ru)- Each line depicts the ratio versus a for a fixed f3. 

We have proved in Proposition 4 that if (oik, is in B 2 f]Ru, then the ratio L k+2 /L\ G [1, 2]. However 
if we analyze each level of fusion, it can be seen that the total error probability decreases exponentially 
fast from B 2 fl Ru to B 1 (Proposition 2). Proposition 6 tells us that the fusion from B 1 to Rc is a bad 
step, which does not contribute significantly in decreasing the total error probability. 

We can now provide bounds for the total error probability at the fusion center. 

Theorem 3: Suppose that (a , /3 ) G R, log N is an odd integer, and there exists k e {0, 1, . . . , logiV — 
1} such that (a fc , @ k ) G B 2 n ifo. 
If is even, then 




2N logL, 



o 1 




) 



^logP^ 1 ^ V2JV login . 



If fc is odd, then 




Proof: If (a m ,/3 m ) E B 2 n R u and m is even, then by Proposition 2, we have 



L. 






2 (fc+l)/2 





where a.j G [1, 2]. 

Let k = log N, we have 
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Therefore, 



\ogV2N 

hgP^ = -[ Yl lo § a * I + v^VlogLo 1 . 

i=i 



Notice that logL 1 > and for each i, < loga^ < 1. Thus, 
Finally, 



logP/^v^logLo 1 . 



log P N X >-\ogV2N + V2N log L- 1 = V2N (log L„ 1 - ^J^ ^J 

If (a m , /3 m ) E B 2 n Ry and m is odd, then by Proposition 6 we have 

L m +2 — a L m +\ 



for some a G [1/2, 1]. 

It follows from Proposition 4 that 



Lk+2 — &kL\ 



for A; = 0, 2, . . . , m — 1, m + 2, . . . , log TV — 2 and some a fc G [1, 2]. Therefore, 



(fc-l)/2 



2 (fe-l)/2 





i=l 



where 1 < a* < 2 for i = 1, 2, . . . , (k - l)/2 and 1/2 < a < 1. Hence, 



f '(logAT-l)/2 



/log x/f \ 



nl 2 (logiV-l)/2 



n 

V 1=1 ) 



and so 



logiV = -loga- 



logo. 



i=l 



Notice that logL 1 > and for each i, < logaj < 1 and — 1 < log a < 0. Thus, 

logP^< y^logL^ + l. 

Finally, 



log /7V 



log P N ' > - log + J ^ log L " 1 = J | | log L 1 



Finally, by combining all of the analyses above for step-wise reduction of the total error probability, 
we can write general bounds when the initial error probability pair (a , /3 ) lies inside B m , where m/1. 



Theorem 4: Suppose that (a , f3 ) G B m , where m is an integer and m/1. 
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If logiV < m — 1, then (Theorem 1) 



N (logLo 1 - ^) < logP/ < N\ogL-\ 



If log N > m — 1, and log iV — m is odd, then 



y/2 m ~ 1 N logLo 1 - S < logP/ < v / 2^ =I iVlogL- 1 . 



V2 m - 1 N 

If log N > m — 1, and log — m is even, then 



V2 m " 2 A^ log L 



1 log V2 m - 2 iV 
" v / 2 m " 2 iV 



< logP^ 1 < V2™N\ogL \ 



The proof uses similar arguments as that of Theorem 2 and it is provided in Appendix A. 



C. Invariant region in B\ 

Consider the region {(a, (5) G U\(5 < y/a and > 1 — (1 — a) 2 }, which is a subset of B\ (see Fig. [9]). 
Denote the union of this region and its reflection with respect to (3 = a by S. It turns out that S is also 
invariant. 

Proposition 7: If (ak ,(3k ) G S 1 , then (a fc , G S 1 for all /c > /c . 
The proof is given in Appendix B. 



CO. 



1 

0.8 
0.6 
0.4 
0.2 

0. 





—B 1 
---S 


/ '/ X 

r ' x 

' X 







Fig. 9. Invariant region S (between dashed lines) lies inside £?i (between solid lines). 

We have given bounds for P/v, which is (twice) the total error probability. It turns out that for the case 
where (a , (3$) G S, we can bound Type I and Type II errors individually. 
Proposition 8: If (3k) G S, then 

1<^<4 



and 



1<%?<4. 
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Remark: It is easy to see that as long as the system stays inside B\, then in a similar vein, these ratios 
0^+2/°^ an d are lower bounded by 1 and upper bounded by a constant. But recall that B\ is 

not an invariant region. Thus, it is more interesting to consider S. 

Proofs are omitted because they along similar lines to those in the other proofs. As before, these 
inequalities give rise to bounds on sequences {a k } and {(5k}- For example, for {at}, we have the following. 

Corollary 2: If (ato, A)) e ^ men f° r even ^ we nave 

2 fc/2 (W - ^2) < log^ 1 < ^logao 1 . 

We have derived error probability bounds for balanced binary relay trees under several scenarios. In 
the next section, we will use these bounds to study the asymptotic rate of convergence. 



V. Asymptotic Rates 

The asymptotic decay rate of the total error probability with respect to N is considered while the 
performance of the sensors is constant is the first problem to be tackled. Then we allow the sensors to 
be asymptotically crummy, in the sense that ao + /?o — > 1. We prove that the total error probability still 
converges to under certain conditions. 



A. Asymptotic decay rate 

Notice that as iV becomes large, the sequence (3k)} will eventually move into the invariant region 
R at some level and stays inside from that point. Therefore, it suffices to consider the decay rate in the 
invariant region R. Because error probability bounds for trees with odd height differ from those of the 
even-height tree by a constant term, without loss of generality, we will only consider trees with even 
height. 

Proposition 9: If L = a + (3 is fixed, then 

logP/~logLoW 

Proof: If Lq = ao + (3 is fixed, then by Theorem 1 we immediately see that Pjy — > as — > 00 
(logP^ 1 — > 00) and 

1 logy/A < logP.y 1 <x 
logLoW " logLoW " 

In addition, because logy/N /\/N — > 0, we have 

logLoW 

which means 

logP^ ^ log L^VN. 

m 

This implies that the convergence of the total error probability is sub-exponential, more precisely, the 
decay exponent is essentially y/N. 

Given L e (0, 1) and e G (0, 1), suppose that we wish to determine how many sensors we need to 
have so that P N < e. The solution is simply to find an N (e.g., the smallest) satisfying the inequality 

Viv(logL -'- 1 ^)>-log £ . 

The smallest N grows like 6((loge) 2 ) (cf., [|30ll , in which the smallest N has a larger growth rate). 
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B. Crummy sensors 

In this part we allow the total error probability of each sensor, denoted by L^ N \ to depend on iV but 
still to be constant across sensors. 

If is bounded by some constant L e (0, 1) for all N, then clearly Pjy — > 0. It is more interesting 
to consider — > 1, which means that sensors are asymptotically crummy. 

Proposition 10: Suppose that Lq N ^ — > 1, or specifically write = 1 — r] N with 77^ — > 0. Then, 



P N ->■ if and only if rj N = u(l/VN). 
Proof: For sufficiently large N, 

1 /r( Ar )\-l 
We conclude that P N — > if and only if 

^log^* ) -1 ->• 00. 

Therefore, 



iVlog^V^-v/iVMl-^). 

But as rr — > 0, — log(l — rr) ~ x/ln(2). Hence, P N — > if and only if r] N VN — )■ 00, often written 

ri N = u(l/VN). 



Now suppose that ci/yN < r/ N < c 2 /vA^. In this case, for large iV we deduce that 

ci < logP^ 1 < c 2 , 

or equivalently, 

2 _C2 < Ptv < 2" C1 . 

Finally, if 77^ = o(l/v / 7V) (i.e., 77^ converges to strictly faster than l/y/N), then Pjv — > 1- 



VI. Conclusion 

We have studied the detection performance of a balanced binary relay tree of sensors and fusion nodes. 
We precisely describe the evolution of error probabilities in the (a, (5) plane as we move up the tree. 
This allows us to deduce error probability bounds at the fusion center as functions of N under several 
different scenarios. These bounds imply that the total error probability converges to sub-linearly, with 
a decay exponent that is essentially y/~N. In addition, we allow all sensors to be asymptotically crummy, 
in which case we deduce the necessary and sufficient condition for the total error probability to converge 
to 0. All our results apply not only to the fusion center, but also to any other node in the tree network. 
In other words, we can similarly analyze a sub-tree inside the original tree network. 

Needless to say, our conclusions are subject to our particular architecture and assumptions. A series 
of questions follows: Considering the same tree configuration, how does the error probability behave if 
prior probabilities are not equal? Considering balanced binary relay trees with sensor and/or connection 
failures, how would the error probability behave? More generally, what can we say about unbalanced relay 
trees? We also would like to apply our methodology to other configurations, such as a tandem structure. 
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Appendix A 
Proof of Theorem 4 

If log N < m — 1, then this scenario is the same as that of Theorem 1. Therefore, 



TV (logLo 1 - < logP/ < AHogLo 1 . 



If logiV > m — 1 and logiV — m is odd, then it takes (m — 1) steps for the system to move into B 1 . 
After it arrives in B\, there is an even number of levels left because log iV — m is odd. 
By Proposition 2, we have 

r _ ~ 7 

J k 



Lk+i — o-kLl 



for k — 0, 1, . . . , m — 2 and some a k e [1,2], and in consequence of Proposition 4, 

Lk+2 = Q>kL\ 

for k = m — l,m — 3, . . . , log N — 2 and some a k e [1,2]. Thus, 



fc-(m-l) 



'm-l+ 2 



/-/■• = Yl "' L 



m— 1H ^ 



J 



i=l 



where a; G [1, 2]. 

Let k = log iV. Then we obtain 



N 



'log v / 2 m - 1 Af \ / log v / 2 m - 1 Af 

n ^U° logv/2 ^= n 



i=l / \ i=l 



and so 



/ log y/2 m ~ 1 N \ 

logP/ = - I £ log a, j + V2— liVlogLp 1 . 
Note that logLp 1 > 0, and for each i, < log a; < 1. Thus, 



logP^ 1 < v / 2^ =I 7VlogL- 1 . 

Finally, 

(log V2 m - 1 N\ 
I +v / 2^iVlogL 1 

= - log \ / 2 m ~ 1 N + v / 2 m " 1 iV log Lq 1 

= v^jvfio g ^- io g 7 gy 

For the case where log iV — m is even, the proof is similar and it is omitted. 
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Appendix B 
Proof of Proposition 7 

Without loss of generality, we consider the upper half of S, denoted by Su- As we shall see, the image 
of Su is exactly the reflection of Su with respect to the line /3 — a (denoted by Sc). We know that 
S u := {(a, (3) e U\/3 < \/ol and /5 > 1 — (1 — a) 2 }. 

The image of Su under / can be calculated by 

(a',/3 / ) = /(«,/3) = (l-(l-«) 2 ,/3 2 ), 
where (a, /3) eU. The above relation is equivalent to 

Therefore, we can calculate images of boundaries for R u under /. 
The image of the upper boundary (5 < ^fa is 



y/p < \/l-Vl-a' 

ol > 1 - (1 - /3') 2 , 
and that of the lower boundary /3 > 1 — (1 — a) 2 is 

v^ 7 > l - (l - (l - VT^))\ 

i.e., 

a' < ^J'. 

The function / is monotone. Hence, images of boundaries of Su are boundaries of Sc- Notice that 
boundaries of Rc are symmetric with those of R u about f3 — a. We conclude that S is an invariant region. 
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